OPTIMAL SYSTEMS OF FUNDAMENTAL ^-UNITS 
FOR LLL-REDUCTION 



L. HAJDU 

Abstract. We show that a particular parameter plays a vital 
role in the resolution of S'-unit equations, at the stage where LLL- 
reduction is applied. We define the notion of optimal system of 
fundamental S-units (with respect to this parameter), and prove 
that such a system exists and can be effectively constructed. Ap- 
plying our results and methods, one can obtain much better bounds 
for the solutions of S'-unit equations after the reduction step, than 
earlier. We briefly also discuss some effects of our results on the 
method of Wildanger and Smart for the resolution of S'-unit equa- 
tions. 



1. Introduction 

S'-unit equations play a central role in the theory of Diophantine 
equations. On the one hand, there are many Diophantine problems 
which naturally give rise to S'-unit equations. On the other hand, 
several types of classical Diophantine equations (such as e.g. norm form 
equations, discriminant form equations, index form equations) can be 
reduced to such equations. Here we only refer to the papers, survey 
papers and books [16], [23], [7], [6], [17], [26], [11] and the references 
given there. 

It is well-known (see e.g. [23], [7], [T7], [8] and the references there) 
that under general conditions the number of solutions of S'-unit equa- 
tions is finite. However, if the number of variables is more than two, 
then no bound is known for the solutions themselves. Moreover, a con- 
jecture of Cerlienco, Mignotte and Piras [4] states that such equations 
are algorithmically unsolvable. In case of S'-unit equations in two vari- 
ables, by the help of Baker's method the solutions can also be bounded 
(see e.g. [15], [23], [7], [17], [18], [19] and the references given there). 
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Further, based on LLL-type results of de Weger (cf. p9j) and an enu- 
meration method of Wildanger f3D|, there exists an efficient algorithm 
for determining all solutions of any particular equations in this case; 
see [25] . 

The practical solution of S'-unit equations in two unknowns (and 
also of some other types, e.g. of index form equations) consists of three 
main steps (see e.g. [21], [13], [H], [30], [25], [10], [1], and the references 
given there). First, applying Baker's method an initial upper bound 
Cini is obtained for the unknowns (which are in the exponents). For 
the best known bounds see [20] in the complex and [21] in the p-adic 
case, respectively, and also the references given there. Roughly speak- 
ing, at this principal stage the "large" solutions are excluded. Though 
from the theoretical point of view this first step is the deepest and 
most important one, for practical purposes (i.e. for listing explicitly all 
solutions) the Baker-type bound is (typically) too large. Hence before 
going for the solutions, some further reduction is needed. In the sec- 
ond step applying some variants of the LLL-algorithm, one can reduce 
the bound Cini considerably, to get a much smaller bound Cred for the 
solutions (see e.g. [2S] and the references given there). This stage is 
frequently referred to as getting rid of the "medium" size solutions. 
However, if the number of variables (the rank of the S'-unit group) is 
"large", even this reduced bound can be too high for getting all solu- 
tions, at least if one would like to simply apply some primitive listing 
algorithm. So to find the solutions effectively, some clever enumera- 
tion is needed. Such an algorithm has been worked out by Wildanger 
[5D] in the complex case, and extended and adopted by Smart [25] for 
the S'-unit case. The application of these algorithms as a third step 
(ideally) leads to the explicit solution of the original S'-unit equation. 

In this paper we focus on the second stage of this procedure. It 
turns out (in fact it is widely known already) that in pushing down the 
reduced bound Cred as much as possible, a certain parameter plays an 
important role. The aim of the present paper is to show that there is 
an optimal choice for this crucial parameter, and to give a method to 
actually find it. Hence ultimately we are able to get much better bounds 
Cred in case of particular equations than earlier. Having established 
our results, we briefiy discuss their possible effects on the method of 
Wildanger [30] and Smart [2S] for the practical resolution of S-unit 
equations. We note that the results presented in the paper lead to 
certain improvements of the method. 

To formulate our results clearly, we need some preparation. For this 
purpose we introduce some (standard) notation. 
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1.1. Valuations. Let K be an algebraic number field of degree n, and 
let Mr be the set of places on K. Choose from every place v G a 
valuation in the following way. If v is infinite and corresponds to 
an embedding o" : K — C then for every a G IK let 




, if cr is real, 



2 



if a is complex. 



Further, if v is finite and corresponds to a prime ideal P of K, then for 
every a G K put 



Jo, if a = 0, 

1 Norm(P)~°''^^("), otherwise. 



By these choices we have the product formula, that is for every a G K, 

(1) H \a\, = l 

holds. Let 5* = {vi, . . . , f^} be a finite subset of Mr, containing all the 
infinite places, and write 

Us = {e: \e\y = 1 for all v e Mk\S} 

for the set of S-units. As is well-known. Us is a finitely generated group 
of rank s — 1, containing the unit group of the ring of integers of K. 

1.2. S-unit equations. Let ai,a2 be fixed non-zero elements of K. 
Consider the so-called S'-unit equation 

(2) aiXi + a2X2 = 1 

in two unknowns Xi,X2 G Us- It is well-known (see e.g. [7J) that if 
Xi, X2 is a solution to equation ([2D then we have 

(3) X, = el'''e\''' . . . e^iV^ for ^ = 1, 2. 

Here ei, . . . is a fundamental system of S'-units, and is a root 
of unity in K. Put Bi = max for i = 1,2 and B = max(i?i, B2). 

l<j<s—l 

Clearly, without loss of generality we may assume that B = Bi. From 
(I3D, for all f G Mk 

s-l 

\og\xi\y = y^&ijiog|g-,i„ 
3=1 
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(7) |1 - a2X2\vt < |ai|t,texp 



holds. In particular, we have 

/logjeil^i ••• log|£s-l|^,A / 61,1 

\log|£i|^^ ... log|e,_i|^^/ 

From the above equality (jl]) one easily gets that 

(5) i?i < C* max I log |a;iL- 1 

l<j<s ' ' ^ ' 

with some constant C*, depending only on our fundamental system of 
S'-units. (The constant C* is the crucial object from the point of view of 
the present paper; we shall come back to this point a bit later.) Then, 
by the help of the product formula ([T]), using standard arguments (see 
e.g. [Z], [H]) we immediately get that 

(6) 

holds for some t G {l,...,s}. Hence from ([2]) we get 

[s - l)C* 

Now using a Baker-type result (e.g. Matveev i20j if vt is infinite and 
Yu [31j if Vt is finite, respectively) we get something like 

(8) |1 - a2X2\vt > exp(-Colog(52)). 

Here Cq is some constant depending only on ai,a2, K, S. The inequali- 
ties ([7]) and ([8]) by our assumption B = Bi yield an initial upper bound 
Cini for B. 

1.3. The importance of C* . To get the reduced bound C^ed by the 
LLL-algorithm, one starts from inequality ([7]) (or a variant of it) to- 
gether with the already known information B < Cini- Interestingly, 
the finally obtained reduced bound Cred depends very heavily on the 
constant C*; the dependence is closely linear. This phenomenon should 
hopefully be (heuristically) clear from the examples presented in this 
paper, but already is well-known for experts for a long time; see e.g. 
the remarks of Tzanakis and de Weger [28] pp. 239-240, and Smart 
[21] p. 823. Hence it seems to be worth to try to keep C* as small 
as possible. Apparently, so far this point remained more or less un- 
touched, and the calculation of C* is usually done in a rather casual 
way. Namely, in all occurrences in the literature the standard choice is 
to take something like 

(9) C* := min Hi^T^II 
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(with certain refinements at some places). Here Rj is the (invertible) 
matrix obtained by deleting the j-th row of the matrix at the left hand 
side of dl]), and stands for the row norm of a A;i x k2 type real 

matrix A = {aij)i<i<ki, i.e. \\A\ \ = max Yl ^jjl- Using e.g. Cramer's 

i<j<fe2 l<i<A;i j—i 

rule, one can easily see that this choice of C* is appropriate to have 

iD- 

However, as it turns out, this choice of C* can be rather far from 
being optimal, which results in a much worse value for Cred than possi- 
ble. In the second section we show a very simple way to get an instant 
improvement upon the above choice. Further, we show that in fact the 
best C* value exists, and depends only on the choice of the fundamental 
system of 5'-units. Finally, we prove that one can explicitly determine 
an optimal fundamental system of S'-units (yielding the best value for 
C*), and we give a (relatively) efficient algorithm for finding such a 
system. As the steps of our arguments and methods are connected in 
a rather organic way, we do not start with listing the main theorems, 
we prefer to formulate our results in a "linear" way. However, in order 
not to break the presentation, we give the proofs in a later section. In 
Section [3] we give an algorithm which finds an optimal system of funda- 
mental S'-units for LLL-reduction. In the fourth section of the paper we 
provide some numerical examples, including the bounds Cred obtained 
by the old method and by the new one. In the fifth section we give the 
proofs of our results. The sixth section is devoted to a brief discussion 
about the effects of our results on the method of Wildanger [30] and 
Smart [25] for the practical solution of S'-unit equations. Finally, in the 
Appendix on the one hand we outline the reduction methods, and on 
the other hand we indicate how one can adjust the method developed 
in the paper if the valuations are not chosen in the "standard" way. 

Note that our method can be adopted to the case where in ([2]) not 
a full system of fundamental S'-units are involved, or the S'-units form 
only an independent system. 

Finally, we mention that a similar type investigation has been per- 
formed about Mordell-Weil bases of elliptic curves by Stroeker and 
Tzanakis [27], to reduce the final bound for the integral solutions of 
elliptic equations. 

2. Optimizing C* 

We keep our notation from the previous section. Further, let F = 
(ei, . . . , Es-i) be a system of fundamental S'-units of K, and define the 
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s X [s — 1) matrix Rp by 



/log 


l^^l |di 


. log|es_i|j,A 


Vlog 


\^l\vs ■■ 


. log\es-i\vJ 



Let R'p be any (s — 1) x s matrix such that 
(10) R'p ■ Rp = E(^s-i)x{s-i), 

with the identity matrix of size (s — 1) x (s — 1) on the right hand 
side. The importance of the matrices R'p becomes clear in view of the 
following simple observation: starting from (jlj), by ( fTOl) we get 

bi,i \ /log|xi| 

A,s-i/ \logFil 

Hence in ([5]) we can take C* to be ||-R^||, with any matrix R'p. So we 
define the norm N{F) of the system F by 



N{F) := min||i?^ 



R'p 



F\ 



where R'p runs through the matrices for which (|TOl) is valid. (We shall 
see later that the minimum does exist.) The system F is called optimal 
for LLL-reduction, if N{F) is minimal among all choices of fundamental 
systems of S-units. As we shall also see, the minimum of N{F) also 
exists. Hence we put 

(11) C* := minA^(F), 

where F runs through all systems of fundamental S'-units. Then we 
have ([5]), of course using the optimal system F in Q. In order to 
compare this choice of C* with the earlier one, we write Noid{F) for 
the choice of C* in ([9]). Note that ([9]) just means that we unnecessarily 
restrict ourselves to matrices R'p having a constant zero column. 

In the following proposition we describe the structure of the matrices 
R'p,, for fixed F. 

Proposition 2.1. Let R'p be a matrix for which ( ITOi) is valid. Then 
for each i G {1, . . . , s — 1} there exists a G M such that the i-th row 
of R'p is of the form Wj — -Uj ■ 1 with Wj = (wj^i, . . . , Wt^s-i, 0), where 1 
is the s-tuple with all entries equal to 1, and 

••• \ /logjeil^, ... logle^^il^, ^ 

... Ws^i,s-i/ \log|ei|„,_i ••• logks-i|i,,_ 
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The next result provides a simple tool to calculate N{F) for a fixed 
system F of fundamental S'-units. Before its formulation, we need a 
new notation. Define the central norm |x|c of x G M", x = (xi, . . . , x„) 
in the following way. Let yi, . . . ,yn be a rearrangement of Xi, . . . , x„ 
such that Vi < ■ ■ ■ < Vn- Then put 

n 

(12) \^\c = ^\yi-y^\ 

1=1 

where / = [{n + l)/2j. The number yi is called the center of x. 
Proposition 2.2. Using the previous notation, we have 

N(F) = max Iwjlc- 

l<i<s-l 

Note that the above result already gives a tool to improve upon 
the earlier choice of C* in ([9]). Indeed, if one does not interested in 
finding the optimal system of fundamental S'-units for LLL-reduction, 
but insists on his favorite system F, taking C* to be N{F) rather than 
Noid{F) in ([5]), is already an improvement. By Proposition 12.21 the 
calculation of N{F) takes only a fraction of time. 

The next result (together with its proof) shows that there exists an 
optimal system of fundamental S-units indeed, and such a system can 
be effectively determined. 

Theorem 2.1. For any positive real constant c there are only finitely 
many systems F of fundamental S -units such that N{F) < c. Further, 
all such systems can be effectively determined. 

3. An ALGORITHM TO DETERMINE AN OPTIMAL SYSTEM F 

In this section we outline an algorithm to determine an optimal sys- 
tem of fundamental S'-units for LLL-reduction. The algorithm consists 
of two parts. First, by a heuristic method we (hopefully) obtain the 
best system F, then we check that our choice is best possible indeed. 

In the first step we start from an arbitrary system Fo of fundamental 
S-units. (Such a system can be found e.g. by Magma [3], but one can 
also use KASH [5] or PARI/GP [2]). Then we calculate the values 
of and choose Ui to be the center of Wj (1 < z,j < s — 1) in 
Proposition 12.11 Hence for the corresponding matrix R'p^ we have that 
I iR'pg 1 1 is minimal (see the proof of Proposition 12. 2p . To find an optimal 
system F, in fact we need to find an unimodular (s — 1) x (s — 1) type 
matrix Aq such that N{FoAo) is minimal. 

For any row vector b = (6i, . . . ,bn) write b* = (6i — -u, . . . , 6„ — m), 
where u is the center of b. We call b* the centralization of b. If i? is a 
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matrix then let B* denote the matrix obtained from B by replacing all 
its rows by their centralizations. Observe that for any matrices A, B 
(of appropriate sizes) we have {A - B)* = {A - B*)*. Hence our problem 
reduces to finding the following minimum: 

(13) min||(/li?^J*|| 

where A runs through all the unimodular matrices of type (s — 1) x 
(s — 1). If this minimum is taken at some matrix A', we simply have 
Ao = {A')-'. 

We determine the above minimum with a heuristic algorithm, which 
seems to work well. The algorithm produces a sequence of systems 
Fo, Fi, F2, . . . with A^(Fo) > A^(Fi) > N{F2) > . . . , and terminates at 
some point. Suppose that we have already made i steps, and currently 
we are working with Fj. Let 



/ ^1,1 ■ 




(i) 





belong to the optimal choice, i.e. ||-Rir-|| = N(Fi). Write wf = 
(wf . . . , wl^^g) for the t-th row of R'p {t = — 1). Suppose 

that ll-R^II = X] I'^jil (that is the j-th row of R'p_ yields the row norm 

1=1 

\\R'f. II). To improve upon | |-R^. 1 1 we need to find an unimodular matrix 
A such that for some row a = (ai, . . . , a^-i) of A we have aj 7^ 0, and 
further 

s-l 

i=l 

Indeed, otherwise we cannot "replace" the vector w^*'' by any "shorter" 
one, and hence we cannot improve upon ||i?^. ||. In the first part of the 
algorithm we try to find such a row vector a of a simple shape, and then 
we iterate the procedure. More precisely, we consider the minimum 
f[T^ . however, only for A running through the unimodular matrices 
which are different from the identity matrix only in their j-th row; 
namely, the (j, j)-th entry of A equals 1, and all the other entries in its 
j-th row may assume the values —1, 0, 1 only. Having the minimum at 
A', we define Fj+i by Fj+i = Fi{A')^^, and then repeat the procedure. 
By Theorem 12. II we know that this algorithm terminates, and produces 
a system Fk as output. 

Now we should check that the final system is optimal. (Note that 
this was the case for every example we tested the algorithm for, so 



< II^F. 

I I r-i 
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we think that this is the typical phenomenon.) For this we need some 
preparation; in fact we need to estabhsh two simple properties of the 
central norm. 

Lemma 3.1. The central norm is homogeneous, that is for any x G M" 

and t G M we have \tx.\c = \t\\x.\c- 

Lemma 3.2. Let Ci G M" denote the vector with all coordinates 0, 
except for the i-th one which is 1 {i = 1, . . . ,n — 1) , and put 

T := {ej : z = 1, . . . , n - 1} U {-e^ : i = 1, . . . , n - 1} U {eo, -eo} 

where eg = ei + ■ • • + e„_i. Further, set 

:= {x G M" ■.Xn = and |x|c; < 1} 

where x„ is the last entry of ^. Then H is the convex hull ofT. 

As a corollary of the above two lemmas we get the following state- 
ment, which will be needed in the last step of our algorithm. 

Lemma 3.3. Let F he a system of fundamental S -units, and choose 
an R'p as before. Then for any (integral) unimodular matrix A of type 
(s — 1) X (s — 1) with ||(y4_R^)*|| < N{F) we have that the row vectors 
of A belong to the convex hull of the set 

{±iV(F)bi,...,±iV(F)bJ, 

where bj is the i-th row of Rp {i = 1, ■ . ■ , s) . 

Now we can outline the final part of our algorithm. Having the 
output Fk by the first stage of the procedure, we calculate N[Fk) (by 
the help of Proposition 12.21) . Then using Lemma 13.31 we get an up- 
per bound Co for each entry of the possible unimodular matrices A. 
Typically, this bound is very small (around at most 2-3). We keep the 
notation introduced above Lemma [3. II In particular, we assume again 
that the value of | |-R^ 1 1 belongs to the j-th row of R'p . Then we check 
all row vectors a = (ai, . . . , ag-i) such that aj > and all the entries 
of a are at most cq in absolute value. If for all such vectors a we have 
l^-^Klc — ^(-^'s)' then the value of N(Fk) cannot be improved and 
Fk is an optimal system. Otherwise, for every appropriate row vectors 
a we check all unimodular matrices with entries at most cq in absolute 
value, containing a as a row. In case we can improve upon N{Fk), 
we switch back to the first part of the algorithm (with the improved 
system), then return to this second part later on, and so on. By The- 
orem 12.11 we eventually get an optimal system of fundamental S'-units. 
We mention that the second part is much more time consuming than 
the first one, at least if s is "large". (It is not surprising in the light 
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of the observation that the number of row vectors a to be checked is 
Co(2co + 1)'^"^ at this stage.) However, note that in all cases we en- 
countered the second part has never been necessary in the sense that 
the system Fj, obtained by the first part of the algorithm has already 
been optimal. 

We conclude this section by two remarks. First we note that our 
simple heuristic algorithm works so well is very probably due to the 
fact that the systems of fundamental units provided e.g. by Magma 
are already LLL- reduced, hence the initial system is "closely" optimal. 
The other thing we mention is that our algorithm is not optimized, it 
probably can be improved. Further, maybe there are other (possibly 
more efficient) ways to get the optimal system. For example, after hav- 
ing an original system Fq one can search for the minimizing unimodular 
matrix A in f|T3l) in the following way. Let C denote the lattice spanned 
by the vectors w^'''' (t = 1, . . . , s — 1). If Fq is not optimal then there is 
some unimodular matrix A such that for all row vectors a of A we have 

s 

'Yj Wt\ < N{Fo) where := {a[, . . . , a'g) = (ai?^^)*. This implies that 

t=i 

Wa is a vector of the lattice C of Euclidean length less than y/sN{FQ). 
Hence all the appropriate row vectors a can be efficiently determined 
by the algorithm of Fincke and Pohst j9j. Having the set of possible 
vectors a, we can build up the appropriate unimodular matrices A, and 
we can find the minimum 0131) . However, as we mentioned, the original 
algorithm was efficient enough for our present purposes, so here we do 
not take up the problem of optimizing the method. 

4. Examples 

In this section we present some examples to illustrate our method 
and also to give some numerical evidence why the parameter C* is 
so important. In fact we have worked out several examples, but we 
give here only four, out of which three can be found in the literature. 
We exhibit a "proper" S'-unit equation (i.e. with a finite valuation 
involved) and three pure unit equations. Among the latter ones the 
first two correspond to totally real fields, while the last one to a totally 
complex field. These examples are of larger size. We mention that in 
case of smaller examples, our algorithm worked with the same efficiency, 
but much faster. However, naturally in the "small" cases the gain using 
the new values of C* is certainly smaller. 

In each example we consider an S-unit equation of the shape 



(14) 
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in xi,X2 G Us, for some particular choice of K and S. Just as in ([3]), 
we write 

(15) Xi = el''°e\''' . . . e'^f ^ for z = 1, 2. 

Note that here the system ei, . . . ,es-i is certainly not fixed during 
the examples, in fact we use two different systems in each example 
(corresponding to different values of C*). 

In every example, starting from an initial system Fq of fundamen- 
tal S'-units, using our algorithm we determine an optimal system F^. 
(Note that at the examples from the literature we use the system Fo 
from the corresponding papers, otherwise we used Magma [3] to get an 
Fq.) Remark that the first stage of our algorithm typically takes a few 
seconds (up to a few minutes when s is larger) to terminate and provide 
an optimal system F^. However, to check that Fk is optimal indeed by 
the second part of the algorithm takes much more time (around 30 
minutes for "large" s). An optimized and more sophisticated version 
of the algorithm would probably work better, but we do not take up 
this question here. 

Description of the tables. We provide a table for each example, 
containing several data. We give the values of the reduced bounds 
Cred corresponding to the choices C* = Noid{Fo), N{Fo) and N{Fk), 
respectively. Naturally, in the first two cases we use the system Fq in 
(fT5|) . while in case of C* = N{Fk) the system F^ is used. To execute 
the reduction steps, we used Lemmas 17.11 and 17.21 given in the Appen- 
dix. Note that these reduction lemmas have very many variants in the 
literature. We chose these ones because they are relatively simply for- 
mulated, and they are appropriate for our present purposes. To make 
the reduction, when there was no available initial upper bound for B, 
as it is not an important point from the point of view of the present 
paper, instead of going through Baker's method we just started with 
the (plausible) bound B < 10000 =: Cini- 

In the tables we provide the ratios of the C* values, as well (more 
precisely the ratios corresponding to the actual C* and Noid{Fo)). We 
also indicate the ratios of the Cred values (that is, the ratios corre- 
sponding to the actual Cred and the reduced bound corresponding to 
the choice C* = NoidiFo)). Note that these ratios are remarkably close 
to each other, which shows that the dependence of Cred is very nearly 
linear in C* . This phenomenon is not surprising in view of Lemmas 17.11 
and 17.21 the new lower bound for B in each iteration is almost linear 
in 1/Ci in both cases - and 1/Ci is linear in C*. (For the details cf. 
subsection 17. II of the Appendix.) 
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Finally, we introduce an indicator called "domain ratio" to compare 
the remaining domains to be checked after the reduction. This indica- 
tor is defined in the natural way, i.e. as 

2s-2 

J -L 1 \ 
red 



Xi':' + 1 



Here the constants Cred are the reduced constants (C^g^ corresponds 
to the choice C* = Noid{Fo) and C^^^ to the actual choice), and the 
exponent 2s — 2 by (ITS!) is just the number of variables (as bi^ and 62,0 
do not really count). Note that in our examples this indicator shows 
that the size of the domain to be checked using the new method is a 
tiny fraction of the size of the domain remaining by the old method. 

Example 1. This example is from p5]. Let K = Q{'&) where + 1 = 
0. Let S be the set containing the infinite valuations of K, and a finite 
valuation corresponding to the prime ideal P = (it) with vr = 1 — -(9. 
We have N^/Q{n) = 2, and a fundamental system of S'-units is given 
by 

e^ = ^^ + ^^ + 52 = -^92 _ # _ 53 = 1 + # - 54 = vr = 1 - ^9 

(see [25]). Note that the element Eq = — -(9^ generates the sixteen roots 
of unity of K. Put 

A simple calculation shows that Noid{Fo) = 1.442695 .... Further, we 
also have A^(-Fo) = 1-442695 .... Then, executing the above algorithm 
starting with Fq, in three steps we get a new system of fundamental 
S'-units F3 given by the transformation 

iooo\ /iooo\ /lo o\ /lo 

1 1 / 1 1 I 1 1 -1 -1 )_/ 1 1 -1 -1 

0010) loiio) loo 1 o)~\ii 0-1 

0111/ Voiii/ Voo 1/ V22-1-1 

Actually, we have 

A simple calculation gives that N{Fs) = 0.931871 .... The second part 
of our algorithm reveals that -F3 is an optimal system of fundamental 
S- units. 

Then, using Lemmas 17.11 and 17.21 from the Appendix, we get the 
reduced bounds in the table below. Note that we need to perform the 
reduction for each choice of the valuations, and that the worst case 
belongs to the finite valuation in S. We also mention that working 
with Fq, we used the initial upper bound B < 1066, from [25j. In 
case of F3, using the inverse of the basis transformation matrix, from 
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this we have B < 3198. We started our reduction with these bounds. 
We summarize the resuhs of our computations in Table 1. For the 
definitions of the entries of the table (and also in case of the other 
tables) see the preceding subsection. 





using Fq and Noid{Fo) 


using Fq and A^(-Fo) 


using F3 and N{F3) 


c* 


1.442695... 


1.442695... 


0.931871... 


Cred 


1031 


1031 


651 


C* ratio 


1 


1 


0.645923... 


Cred ratio 


1 


1 


0.631425... 


domain ratio 


1 


1 


0.025325... 



Table 1 



Example 2. The data of this example is from [TU] (the authors con- 
sidered a different equation). Let K = Q{'&) where — + •t}'^ + 
66^6 + ^5 _ 9g^4 _ 7^3 ^ 37^2 + 12^9 + 1 = 0. Let S be the set of 

infinite valuations of K. An integral basis of K is given by tui, . . . , ujiq 
with LUi = -d^'^ (i = 1, . . . , 9) and 

9 + 27^9 + 43^?2 + 20^93 + 37^94 + 5^?5 + 32^96 + 3^7 ^ 26^9^ + d'^ 
^10 = ^7 . 

The coordinates of a fundamental system ei,...,eg of S-units with 
respect to this integral basis is given by 

[21, 107, 192, -5, -120, -40, 84, 20, 30, -60], 

[16, 99, 139, -56, -113, -7, 56, 9, 14, -30], 

[10, 4, 65, 197, 85, -110, 56, 34, 50, -90], 

[21, 35, 196, 346, 94, -206, 129, 66, 97, -177], 

[0, -53, -31, 200, 145, -90, 14, 24, 35, -60], 

[8, 24, 40, 33, -1, -27, 25, 10, 15, -28], 

[15, 13, 118, 248, 78, -143, 84, 45, 66, -120], 

[0,1,0,0,0,0,0,0,0,0], 

[4, 19, 42, 0, -26, -8, 17, 4, 6, -12], 
respectively (see [IQ]). Put 

-^0 = (^1) ^2, £^3; ^4, £^5) ^6) ^7) ^8) ^g)- 
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We have Noid{Fo) = 2.285921... and A^(Fo) = 1.564168.... Then 
by the above algorithm starting with Fq, in three steps we get a new 
system of fundamental S'-units given by the transformation 



10 








10 








10 








10 1 


-1 


-1 


10 








1 











1 








1 











1 



lOOOlOOOOx 

loooooooN 



/o 

001000000 
000100000 
000010000 
000001000 
000000100 
VOOOOOOOIO/ 
NOOOOOOOOl/ 








-1 




1 
1 




lOOOOOOOOx 

oloooooooN 

001000001 
000100000 
000010000 
000001000 
000000100 
000000010 / 
.000000001/ 



/ 1 1 

/ 10 

10 

10 1 

10 

1 



\ 

\0 



\ 

1 

-1 




?/ 



Actually, we have 



Fs = {e, , 62, 63, £4, EiE^, E^Eq, £4 Eq, Ej, E^E^ 

By a simple calculation we get that NlF^) = 1.209236 .... Using the 
second part of our algorithm we obtain that F3 is an optimal system of 
fundamental S-units. In this case we performed the reduction starting 
with the initial bound B < 10000 =: Cmi both with Fq and with F3. 
The results of our calculations are summarized in Table 2. 





using Fq and Noid{Fo) 


using Fq and N{Fq) 


using F3 and NlF^) 


c* 


2.285921... 


1.564168... 


1.209236... 


Cred 


2079 


1416 


1011 


C* ratio 


1 


0.684261... 


0.528993... 


Cred ratio 


1 


0.681096... 


0.486291... 


domain ratio 


1 


0.000996... 


0.000002... 



Table 2 



Example 3. This example is from [30j. Let IK = Q(i9) where 

1}^ + - 8^9^ - 7# + 2W^ + 15^9^ - 20^^ - 10^9^ + 5^9 + 1 = 0. 

Note that IK is the maximal real subfield of the cyclotomic field Q(Ci9)- 
Let S be the set of infinite valuations of IK. An integral basis of K 
is given by . . . . The coordinates of a fundamental system 
El, . . . ,Es of 5- units with respect to this integral basis is given by 

[1, -4, -10, 10, 15, -6, -7, 1, 1], [0, 3, 0, -1, 0, 0, 0, 0, 0], 

[1, -2, -3, 1, 1, 0, 0, 0, 0], [2, 0, -9, 0, 6, 0, -1, 0, 0], 

[0,1, 0,0, 0,0, 0,0,0], [2, 0,-1, 0,0, 0,0, 0,0], 

[2, 0, -4, 0, 1, 0, 0, 0, 0], [0, -5, 5, 10, -5, -6, 1, 1, 0], 
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respectively; see 
is a typo.) Put 



30] . (Note that in [30j 62 = Si} — is written, which 



-^0 — (^1) ^2, £^3; ^4, £^5; ^6) £^7) ^s)- 



A simple calculation shows that Noid{Fo) = 2.561675 .... Further, we 
also have N{Fo) = 1.872827. . . . Then our algorithm in six steps yields 
a new system of fundamental S-units Fq given by the transformation 



/10 -1 00000\ 

101 ooooooN 












100000 
010000 
001000 
000100 
000010/ 
001/ 



10000000\ 

looooooX 



/o 

00100001 
00010000 
00001000 
00000100 







vooooooio/ 
\oooooooi/ 



10000000\ 
OlOOOOOOX 
00100000 
00010000 
00001100 
00000100 
00000010 / 
1 / 



/ 1 

/ 1 

1 









0\ 

\ 



1 
10 
10 



•0 



\000 00010/ 

\ooo -10001/ 



/10000000\ 
/ OlOOOOOOX 
00100000 
00010000 
00001000 
00000110 
\00000010 / 
\0000000l/ 



That is, we have 



lOOOOOOOx 
OlOOOOOOX 
00100000 
-10111111 
00001000 
00000100 
00000010 / 
00000001/ 





1 

1 
-1 






1 



1 0\ 

\ 

-1 

1 1 

1 
1 
1 , 

-10/ 



^8)^2)^4^8 5 ^1^3 ^4^8 > ^1^3 ^4^5^8 ' 

with N{Fq) = 1.343979 .... The second part of our algorithm verifies 
that Fg is optimal. We started the reduction with the initial bound 
B < 2076 for Fq (see [SU]), while using the inverse of the basis reduction 
matrix we could start with B < 4152 in case of Fg. The data and the 
information derived from them are given in Table 3. 





using Fq and Noid{Fo) 


using Fq and N{Fq) 


using Fq and N{Fq) 


c* 


2.561675... 


1.872827... 


1.343979... 


Cred 


1664 


1210 


824 


C* ratio 


1 


0.731094... 


0.524648... 


Cred ratio 


1 


0.727163... 


0.495192... 


domain ratio 


1 


0.006122... 


0.000013... 



Table 3 



Example 4. This example is new. Let K = Q{i)) where 

^18 + + ^16 + ^15 ^ ^14 ^ ^13 ^ ^12 ^ ^11 ^ ^10 ^ ^9^ 
+ + # + + # + + + ^9 + 1 = 0, 

SO K is the cyclotomic field Q(Ci9) (with ^ig = ■(?). Let S be the set of 
infinite valuations of K. An integral basis of K is given hj , d^"^ . 
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The coordinates of a fundamental system ei,...,es of S'-units with 
respect to this integral basis (obtained by Magma [3j) is given by 

[1,1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,1] 
[0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0] 
[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,1,1] 
[1,1,1,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1] 
[1,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] 
[1,1,1,1,1,1,1,1,1,1,1,1,1,0,1,1,1,1] 
[1,1,1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1] 
[0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0] 
respectively. Put 

-Po = i^l, ^2, £3, ^A-, £^5, ^li ^%)- 

A simple calculation shows that Noid{Fo) = 1.280834. . . and N{Fo) = 
0.936410 .... Then the above algorithm in seven steps provides a new 
system of fundamental S'-units Fj given by the transformation 



10000000\ 
OlOOOOOOX 
00100001 
00010000 
00001000 
00000100 
00000010 / 

1/ 



■0 



/lO OOOOOOx 

/Ol oooooox 

00 100000 

00 010000 

00 001000 

00 000100 

\00 000010/ 

\oo -100011/ 



10000000\ 
OlOOOOOOX 
00100000 
00010000 
00001000 
00000100 
00010010/ 

1/ 



/lO OOOOOOx 
/ 01-11001l\ 














\ 

\o 



10 

10 

10 

10 

1 oy 







10000000\ 

OlOOOOOOX 

00100000 

00011000 

00001000 

00000100 

00000010 / 

00000001/ 



10000000\ 
OlOOOOOOX 
00100000 
00010000 
00001000 
10000100 
00000010 / 
1 / 







/I 






1-11 
1-11 





oil 
oil 



11-1-1 



1-11-111-1 



1 1 
1-11 oil 
10 oil 










10 

'111 

1 









Actually, we have 



-^7 — (£^1^2£^3£^4£^5^6^7^8; ^3 ^4 ^5 ^7 ^^3^i^5^7,^5 



-1^-1 



Sa Err, S 



with N{F7) = 0.67198843 .... By the second part of our algorithm we 
get that F7 is an optimal system of fundamental S'-units. To execute 
the reduction, we used the initial bound B < 10000 both with Fq and 
with F7. The outcome of our calculations is contained in Table 4. 
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using Fq and Noid{Fo) 


using Fq and N{Fo) 


using F-j and N{Fi) 


c* 


1.280834... 


0.936410... 


0.671988... 


Cred 


792 


550 


386 


C* ratio 


1 


0.731094... 


0.524649... 


Cred ratio 


1 


0.694444... 


0.487373... 


domain ratio 


1 


0.002938... 


0.000010... 



Table 4 



5. Proofs 

In this section we give the proofs of our results. 
Proof of Proposition \2. 1\ Noting that the rank of Rp is s — 1 and that 

s 

by the product formula ([T]) we have ^log|£:j|t,. = for each i G 

{l,...,s — 1}, the statement follows from the elementary theory of 
systems of linear equations. □ 

Proof of Proposition \2.S[ In view of Proposition 12. H the only freedom 
we have in the choice of R'p is to choose ui, . . . ,Us-i. Take any i G 
{l,...,s — 1}, and let zi, . . . ,Zs be the rearrangement of the entries of 
Wj such that Zi < ■ ■ ■ < Zs- Writing Wi^s = 0, it is obvious that 



mm >^ \ Wi i — M, ! 



is achieved by the choice = zi, where / = [(s + 1)/2J. Hence the 
statement follows. □ 

Proof of Theorem \2.1\ Let T = {rji, . . . ,?7<i_i) be an arbitrary system 
of fundamental S'-units, and suppose that for the system F we have 
N{F) < c. Then we have Rp = RtA with some integral unimodular 
matrix A of size (s — 1) x (s— 1). Obviously, it is sufficient to "bound" A. 
Observe that in view of Proposition [2^ N{T) and N{F) are exclusively 
ruled by 



AT) 



and Wf 



(T) / 



VwS/ 



respectively (with the obvious notation). Since Wp = A ^Wt, Wt and 

( F) 

Wp are bases of the same lattice. As the last entry of each is 
zero, the restriction |w- |c < c yields that the absolute values of the 
components of w|^'' are all at most c. That is, Wp is a basis of a fixed 
lattice, within a bounded region, and our statement follows. □ 
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Proof of Lemma \3.1[ The statement is a trivial consequence of the def- 
inition f[T^ and of the fact that (using the notation above ffT^ ) by 
symmetry 

n 

|x|c = ^\yi' -yi\ 
1=1 

holds, where I' = [{n + 2)/2\. □ □ 

Proof of Lemma \3.2[ First we prove that any convex linear combina- 
tion of the vectors in T belongs to H. For this purpose, choose ar- 
bitrary non-negative real numbers Aq , Aq , Aj,^, A^, . . . , ^n-i^ Ki-i such 

n-l 

that E (^^ + \~) = 1- Put 

i=0 

n-l 
i=0 

Then for each i G {1, . . . ,n — 1} the i-th entry Xi of x is given by 
Xi = Xq — Xq + Xf — Aj~, while x„ = 0. Hence clearly, 

l^lc = |(A]^ - A^;, . . . , A+_i - A^_i, Aq - Aq )|c. 

Thus 

n— 1 n— 1 

|x|c<$^|A+-Ar|<5^(A+ + Ar) = l, 

and the statement follows. 

Let now x G if be arbitrary, and write x = (xi,...,x„). Let 
...,?/„ be a rearrangement of xi, . . . ,Xn such that yi < • ■ ■ < Vn and 
write Xi = yu(i), where u is the underlying permutation of the indices 
1, . . . ,n. (Note that yu{n) = a;„ = 0.) As usual, put / = [{n + 1)/2J, for 
i = 1,. . . , n-l write A+ = max{y^(i)-yi, 0} and A,~ = max{yi-y^i^i), 0}, 
and set Aq = max{yi, 0} and Xq = max{—yi, 0}. Then on the one hand, 
Xi = (y^ii) - yi) + yi imphes 

n — 1 

x = 5^(A+e, + Ar(-e,)). 

i=0 

Further, on the other hand 

n— 1 n 
i=0 1=1 

as X G ii, and the lemma follows. □ 
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Proof of Lemma \3. 3[ Let be the i-th row of A {i = 1, . . . , s — 1). 
Using Proposition 12.11 without loss of generahty we may assume that all 
entries of the last column of R'p are zero. Hence by < N{F) 

Lemmas 13. II and 13.21 imply that auRp belongs to the convex hull of the 
set 

{±iV(F)eo, ±iV(F)ei, . . . , ±iV(F)e,_i} 

(where the Sj are defined in Lemma [3l2l) . However, then using R'p-Rp — 
-E(s_i)x(s-i) we get that belongs to the convex hull of the set 

{±iV(F)bo, ±iV(F)bi, . . . , ±iV(F)b,_i}, 

where bo = bi + ■ ■ ■ + b^^i. Since by the product formula ([1]) we have 
hg = —bo, the lemma follows. □ 

6. Effects on the method of Wildanger and Smart 

In this section we discuss about the effects of the results in the paper 
on the method of Wildanger [30] and Smart [25] for the resolution of 
S'-unit equations. For this purpose, first we briefly and schematically 
outline the main steps of the method. After that we show that our 
results yield certain improvements of the method. At this point the 
author would like to express his deep thanks to Attila Petho for the 
many fruitful and stimulating discussions and advice about the contents 
of this section. 

6.1. The method of Wildanger and Smart. Now we briefly sketch 
the method worked out by Wildanger ^30j and Smart [25] for the res- 
olution of S'-unit equations. We follow the presentation in [25], with 
certain simplifications. First we need to introduce some notation. 
Let IK and S be as before. For G M with K > 1 write 

{{K, S)) = {aeK: 1/K < < K for all v G S}. 

Let C denote the set of solutions of the S-unit equation ([2]), that is 

£ = {(xi, X2) eUs y< Us ■■ aiXi + a2X2 = !}• 

Further, put 

CH, = {{xi,X2)eC:B<Hi} 
where B is defined after ([3]), and set 

ChAK) = {{XUX2) G Ch^ : xi G {{K,S))}. 

Starting from ([2]) after making the standard steps (see section 1.2) we 
arrive at (jlj). Write Si, . . . for the columns of the matrix at the 
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left hand side of (jll), and X for the vector on the right hand side. Then 
03]) can be rewritten as 

(16) hi,i£i + --- + h^s-i£s-i = X. 

As earher, denote by Cred the reduced bound obtained for B (defined 
after ([3])) after executing Baker's method and the LLL-algorithm. Put 

(17) Kq = maxexp(C^ed| log|£:i|„| H h Cred\ log l^s-iM)- 

Then by the help of flTB]) one can easily get that 

(18) C = Ch,{Ko) 

with Hq = Cred (see Lemma 1 of 
Put 

and 



Si = maxmax(|Q;i|„, -^l^) for i = 1,2 



S3 = maxmin(|a;2 



-1| 



V J 



Now let Ki, Ki^i be real numbers such that max(si, S2, S3, (S3— l)/si) < 
Ki+i < Ki and let Hi G Z. Note that we have Ki^i > 1. For v E S 
define the sets 

Ti^^{Hi,Ki,Ki+i) = \{xi,X2) e Cn.iKi) : | - aiXi - 1|^ < — 

[ 1 + SiAj+i 

1 



T2,y{Hi, Ki, Ki+i) = \ {xi,X2) e CuXKi) ■ 



1 

1 



aiXi 



< 



^ 1 + siKi+i 



T3^y{Hi, Ki, KiJ^i) = {{xi,X2) e CnXKi) '■ I - "2X2 - l|i, < 

012X2 e ((1 + siKi, S))}, 

Ti^^{Hi,Ki,Ki+i) = {{xi,X2) G CH,{Ki) : 



02X2 _ ^ 



aiXi 
012X2 



aiXi 

Further, for j = 1,2, 3, 4 let 



e {{l + s,K„S))}. 



Tj{H„Ki,K,+,) = [jT,,,{H„K„Ki+^] 



Then by Lemma 2 of [25] we have the decomposition 



4 



(19) CnAKi) = CH,^AK^+l)\JT,{H„Ki,K,+,) 
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with 

(20) = C*-C+ 

where C* is the crucial constant investigated in the paper (see ([5])) and 

(21) = max (log (^i^^) . log (^i^^) . 'og(AV,)) • 

For more details and explanation see [30] and [25]. Now what happens, 
is that under certain conditions (which usually hold in the first part 

4 

of the algorithm; see Lemma 3 of [25]) the set |J Tj{Hi, Ki, Ki+i) oc- 

curring in (fT9l) proves to be empty. This is not the case in the later 
part of the algorithm. However, at that stage one can rather efficiently 
use the method of Fincke and Pohst [9J to completely enumerate this 
set. Hence starting from i = 0, by the repeated application of (fT9l) 
the whole procedure can be iterated. Finally we are left only with a 
set CHj{Kj) for some small values of Kj and Hj, which can also be 
enumerated without any trouble. Hence in this way we can get all 
solutions of the original 5'-unit equation ([2]). 

6.2. Some improvements. In this subsection we indicate at which 
points our results may improve upon the above described method of 
Wildanger and Smart. 

By our method presented in the paper we are able to minimize the 
value of Cred- Hence the initial value of Kq in fll7l) can be taken much 
smaller than previously. In particular, one may even use his/her orig- 
inal fundamental system Fq of units. Working with the value N{Fq) 
defined in the paper instead of Noid{Fo) used earlier, we get a better 
Cred than previously. As in this case all the other parameters in f|T7j) 
are unchanged, we already obtain some improvement. 

Furthermore, in the definition fl20|) of -ffi+i the constant C* is used. 
As by our results one can take a much smaller value for this parameter 
than earlier, in each iteration of the algorithm we get a smaller value 
for -ffj+i, and hence the procedure can be made to "converge" faster. 
Here the above remark applies again: using the original non-optimized 
system Fq, but working with the new norm N{Fo), one already gets 
some improvement. 

7. Appendix 

This final section has two distinct parts. In the first subsection we 
provide the reduction lemmas we used to obtain the reduced bounds 
Cred for B, both in the finite and in the infinite case. In the second 
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subsection we show how one can adjust the method for other choices 
of the valuations. (We give our motivation for doing so, as well.) 

7.1. Reduction. There are very many variants of reduction lemmas, 
both in the finite and in the infinite case; we chose a lemma of Smart 
[23] in the p-adic case, and a result of Gaal and Pohst [12] in the 
complex case. Note that these lemmas have to be applied for each 
possible choice of the valuation v & S, and the final reduced bound 
will be the maximum of the bounds obtained for each separate v. 

7.1.1. The p-adic case. To execute the reduction, in the p-adic case we 
use a lemma of Smart [24], which is based upon results of de Weger 
[2U] . For its formulation we need some preparation, in which we follow 
the presentation in [24j with slight modifications. (For more details see 



Let P be a prime ideal in K, corresponding to a finite valuation v 
in S. Suppose that P lies above the rational prime p, and suppose 
that ordp(A) > CiB — C2 with some constants Ci > and C2, where 
A = aiXi = 1 — a2X2 in ([2]). Now assuming that B is not too small 
(otherwise we can derive a much better bound for B than with the 
reduction), we can find /ij G K (z = 0, 1, . . . , s' — 1) with s' = s ot s — 1, 
such that 



where ki E Z with \ki\ < B {i = l,...,s'). Note that these /Xj can 
actually be found (see [2^). 

Let Qp and Kp denote the p-adic closure of Q and the P-adic closure 
of K, respectively. Then we can write Kp = Qp{4>) with some 4> G ^p] 
put riQ := [Kp : Qp]. Further, set 



[23].) 



s' 



0:2X2 = ^oY\_ i^i' 



i=l 




i=l 



Then we can write 



"0-1 



i=0 



where 
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with the corresponding Pji G Qp (i = 0, 1, . . . , no — 1). Let A G Qp such 
that 



ordp(A) = min min (ordp(/?j,i)) . 

l<J<s \ 0<2<no— 1 / 

Then (assuming again without loss of generahty that B is "not too 
small") by a simple calculation, including Evertse's trick (see [28] and 
[21]) we get that 

ordp(A,/A) >CiB-Cs (z = 0, 1, . . . , no - 1). 

Here C3 is a constant which can be explicitly given in terms of C2, A, (p. 
Write 



Ai/ A = Ko^i + ^ kjKj^i, Kj^i eZp, < i < uq- I 
i=i 

with the obvious notation. 

For 7 G Zp and m G Z let 7*^"^ denote the unique rational integer 
such that < 7^"-' < — 1 and 7 = 7*-"^ (mod p"). Further, for n G Z 
set 



/ 1 



0\ 



■-1,0 



\^l,no-l 



and 



^s',no-l ^ 



( \ 



(«) 
^0,0 



^(s'+no)x(s'+no) 



G Z' 



s'+no 



Let £ denote the lattice generated by the column vectors of L over Z, 
and set 



^(^,y) 



min I Ixl I if y = 0, 



mm X 



otherwise. 



The following result is Lemma 5 in [23]. Note that the statement is 
in fact a consequence of Lemmas 3.4, 3.5 and 3.6 of |29j . 
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Lemma 7.1. Using the previous notation, suppose that ordp{A) > 
CiB — C2 with B < Xq. Then i{C,y) > y/s'Xo implies that B < 
iu + Cs)/C^. 



To use this lemma, recall that by % 
\aiXi\v < \ai\y exp 



we have 
-B 



for some v & S. Assuming that v is the valuation occurring in the 
above argument, this yields 

p-Ue,ovd,(arxr) ^ |c,^|^exp 



-B 



where Norm(P) 
simple calculation gives 



s-l)C\ 

pfp and Cp is the ramification index of P. Hence a 



ordp(A) > 



B- 



log I ail 



Jpep\ogp{s - l)C* ) yfpCplogp^ 

Now we can apply Lemma 17.11 with the above inequality, using that 
B < Xq holds for some constant Xq. 

7.1.2. The complex case. To execute the reduction in the complex case 
we use the following result, which is an immediate consequence of 
Lemma 1 of Gaal and Pohst 



Lemma 7.2. Let ^1, . . . , be non-zero real numbers, and let xi, . . . ,Xk 

be integers. Put X = max(|xi|, . . . , \xk\). Suppose that 

\xi^i + ■ ■ ■ + Xfc^fcl < C2 exp(— CiX) and X < Xq 

hold with some positive constants Ci,C2 and Xq. Further, let hi be 
the first vector of an LLL-reduced basis of the lattice spanned by the 
columns of the k x [k + 1) type matrix 

. 

1 . 



/ 1 




.. 



\ 





1 



where H is some positive constant. Then 

\hi\ > v/(A; + l)2'=-i-Xo 

implies that 

log H + log C2 - log Xq 



X < 
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Now we briefly explain how to apply this lemma. For this purpose, 
let f G 5* be an infinite valuation such that |1 — a23;2|t) < C*2 exp(— Ci-B) 
(see ([7])). Then using the inequality | logx| < 2|x — 1| which holds for 
|x — 1| < 0.795, we get 

1^2,1 log l^i I,, H hb2,s^i\og\es-i\v\ < 2C2exp{-CiB). 

Then we can apply Lemma 17.21 with the previous inequality, knowing 
that B < Xq for some constant Xq. 

7.2. Adjusting the method for other choices of the valuations. 

In this subsection we indicate how one can adjust our results for other 
choices of the valuations. The motivation is that if S contains only 
the infinite places (i.e. we are interested in pure unit equations) then 
there is an alternative, also natural choice for the valuations: we simply 
take the absolute value of the real conjugates and the absolute value of 
one from each pair of complex conjugates of elements a G K, without 
squaring in the complex cases. To be more general, keep the previ- 
ous notation, and take arbitrary non-zero rational numbers ri, . . . , r^. 
Choose now valuations such that := \a\^_ for all a G K, where 
\a\v. is the previously defined ("standard") valuation corresponding to 
Vi G S {i = 1, . . . ,s). For simplicity, we do not introduce new nota- 
tion but use the previous one, with the convention that everything is 
adopted for this new choice of the valuations. We have the following 
variants of Propositions 12.11 and 12.21 

Proposition 7.1. Let R'p he a matrix for which ffTOj) is valid. Then for 
each i G — 1} there exists a G M such that the i-th row of R'p 

is of the form Wj - ■ (1/ri, . . . , l/r^) with Wj = {wi^i, Wi^s-i, 0), 
where 

( logl^il^' ... log|£:,„i|^' \ ^ 



• • • wi,, 



\logki|<_, ... \og\es-i\v'^^J 



Note that using the assertion ^ ^ log \ei\^i, = (z = 1, . . . , s — 1) 



• 1 J 



the statement can be proved in a similar manner as Proposition 12.11 
We suppress the details. 

To formulate our last statement, write rj = qi/pi with Pi,qi G Z, 
qi > 0, gcd{pi,qi) = 1 and put ti = q/n where q = gcd(gi, . . . , g^). 

Finally, let be the -tuple such that the first entries of 

equal riWi^i, the next |t2| entries of w- equal r2tfi,2, etc. {i = 1, . . . , s), 
with the convention w,- = = 0. 
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Proposition 7.2. Using the above notation, we have 

max \w^\c 
N(F) = '-^^ . 

q 

s s 

Observe that for z = 1, . . . , s we have = - Yl 

Ui\. Using this assertion, Proposition 17.21 can be verified similarly to 
Proposition 12.21 We omit the details once again. 

Finally, we note that as one can easily see, using Prop osit ions 17 . 1 1 and 
17.21 in place of Propositions 12.11 and 12.21 respectively, all the arguments 
of the paper remain valid also under this general choice of the valuations 
- after making the necessary (but rather obvious) alternations. 

8. Acknowledgements 

The author would like to express his special thanks to Attila Petho 
for the several discussions and advice, in particular concerning Section 
El of the paper. 

References 

[1] Y. BiLU, I. Gaal, K. Gyory, Index form equations in sextic fields: a hard 

computation, Acta Arith. 115 (2004), 85-96. 
[2] C. Batut, K. Belabas, D. Bernardi, H. Cohen, M. Olivier, Users 

Guide to PARI/GP. 

[3] W. BOSMA, J. Cannon, and C. Playoust, The Magma algebra system. I. 

The user language, J. Symbolic Comput. 24 (1997), 235-265. 
[4] L. Cerlienco, M. Mignotte, F. Piras, Suites recurrences lineaires: pro- 

prietes algebriques et arithmetiques, Enseign. Math. (2) 33 (1987), 67-108. 
[5] M. Daberkow, C. Fieker, J. Klitners, M. Pohst, K. Roegner, K. 

WiLDANGER, KANT VI J. Symbolic Comp. 24 (1997), 267-283. 
[6] J.-H. EvERTSE, K. Gyory, Decomposable form equations, in: New Advances 

in Transcendence Theory, A. Baker (ed.), Cambridge University Press, 1988, 

175-202. 

[7] J.-H. EvERTSE, K. Gyory, C. Stewart, R. Tijdeman, S-unit equations 

and their applications, in: New Advances in Transcendence Theory, A. Baker 

(ed.), Cambridge University Press, 1988, 110-174. 
[8] J.-H. EvERTSE, H. P. SCHLICKEWEI, W. M. SCHMIDT, Linear equations in 

variables which lie in a multiplicative group, Annals Math. 155 (2002), 807- 

836. 

[9] U. FiNCKE, M. Pohst, Improved methods for calculating vectors of short 
length in a lattice, including a complexity analysis, Math. Comp. 44 (1985), 
463-471. 

[10] I. Gaal, K. Gyory, Index form equations in quintic fields. Acta Arith. 89 
(1999), 379-396. 

[11] I. Gaal, Diophantine Equations and Power Integral Bases, New Computa- 
tional Methods, Birkhauser Boston, 2002. 



OPTIMAL SYSTEMS OF FUNDAMENTAL S-UNITS FOR LLL-REDUCTION 27 



[12] I. Gaal, M. Pohst, On the resolution of relative Thue equations, Math. 

Comp. 71 (2001), 429-440. 
[13] I. Gaal, A. Petho, M. Pohst, On the resolution of index form equations 

in biquadratic number fields I, J. Number Theory 38 (1991), 18-34. 
[14] I. Gaal, A. Petho, M. Pohst, On the resolution of index form equations 

in biquadratic number fields II, J. Number Theory 38 (1991), 35-51. 
[15] K. Gyory, On the number of solutions of linear equations in units of an 

algebraic number field, Comment. Math. Helv. 54 (1979), 583-600. 
[16] K. Gyory, Resultats effectifs sur la representation des entiers par des formes 

decompo sables. Queen's Papers in Pure and Appl. Math. 56 (1980). 
[17] K. Gyory, Some recent applications of S-unit equations, in: Journees 

Arithmetiques de Geneve 1991, D. F. Coray and Y.-F. S. Petermann (eds.), 

Asterisque 209 (1992), 17-38. 
[18] K. Gyory, Solving Diophantine equations by Baker's theory, A panorama of 

number theory or the view from Baker's garden (Ziirich, 1999), Cambridge 

University Press, Cambridge, 2002, 38-72. 
[19] K. Gyory, K. R. Yu, Bounds for the solutions of S-unit equations and de- 
composable form equations, Acta Arith. 123 (2006), 9-41. 
[20] E. M. Matveev, An explicit lower bound for a homogeneous rational linear 

form, in logarithms of algebraic numbers II, Izv. Ross. Nauk. Ser. Mat. 64 

(2000), 1217 1269. 

[21] A. Petho, Computational methods for the resolution of Diophantine equa- 
tions, Number theory (Banff, AB, 1988), 477-492, dc Gruyter, Berhn, 1990. 
[22] A. Petho, Algebraische Algorithmen, Vieweg, 1999. 

[23] T. N. Shorey, R. Tijdeman, Exponential Diophantine Equations, Cam- 
bridge University Press, Cambridge, 1986. 

[24] N. P. Smart, The solution of triangularly connected decomposable form equa- 
tions. Math. Comp. 64 (1995), 819-840. 

[25] N. P. Smart, Determining the small solutions to S-unit equations. Math. 
Comp. 68 (1999), 1687-1699. 

[26] V. G. Sprindzuk, Classical Diophantine Equations, Lecture Notes in Math. 
1559, Springer, 1993. 

[27] R. J. Stroeker, N. Tzanakis, On the Elliptic Logarithm Method for Elliptic 
Diophantine Equations: Reflections and an Improvement, Experimental Math. 
8 (1999), 135 149. 

[28] N. Tzanakis, B. M. M. de Weger, How to explicitly solve a Thue-Mahler 
equation, Compositio Math. 84 (1992), 223-288. 

[29] B. M. de Weger, Algorithms for Diophantine Equations, Centre for Mathe- 
matics and Computer Science Amsterdam, CWI- Tract 65, 1989. 

[30] K. Wildanger, Uber das Losen von Einheiten- und Indexformgleichungen in 
algebraischen Zahlkdrpen, J. Number Theory 82 (2000), 188-224. 

[31] K. R. Yu, p-adic logarithmic forms and group varieties III, Forum Math. 19 
(2007), 187-280. 



28 



L. HAJDU 



L. Hajdu 

Institute of Mathematics, University of Debrecen 

Number Theory Research Group, Hungarian Academy of Sciences 

AND University of Debrecen 

H-4010 Debrecen, P.O. Box 12, Hungary 

E-mail address: hajdul@math.klte.hu 



